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The properties of Bose-Einstein condensed gases can be strongly altered by tuning the 
external magnetic field near a Feshbach resonance. Feshbach resonances affect elastic 
collisions and lead to the observed modification of the scattering length. However, 
as we report here, this is accompanied by a strong increase in the rate of inelastic 
collisions. The observed three-body loss rate in a sodium Bose-Einstein condensate 
increased when the scattering length was tuned to both larger or smaller values than 
the off-resonant value. This observation and the maximum measured increase of 
the loss rate by several orders of magnitude are not accounted for by theoretical 
treatments. The strong losses impose severe limitations for using Feshbach resonances 
to tune the properties of Bose-Einstein condensates. A new Feshbach resonance in 
sodium at 1195 G was observed. 



Most of the properties of Bose-Einstein condensates 
in dilute alkali gases are dominated by two-body col- 
lisions, which can be characterized by the s-wave scat- 
tering length a. The sign and the absolute value of the 
scattering length determine e.g. stability, internal energy, 
formation rate, size, and collective excitations of a con- 
densate. Near a Feshbach resonance the scattering length 
varies dispersively |Q, || covering the whole range of pos- 
itive and negative values. Thus it should be possible to 
study strongly interacting, weakly or non interacting, or 
collapsing condensates (l| all with the same alkali species 
and experimental setup. 

A Feshbach resonance occurs when the energy of a 
molecular (quasi-) bound state is tuned to the energy 
of two colliding atoms by applying an external magnetic 
field. Such resonances have been observed in a Bose- 
Einstein condensate of Na(i 7 ' = 1, mp = +1) atoms at 853 
G and 907 G |0, 13], and in two experiments with cold 
clouds of 85 Rb(F = 2,TO f = -2) atoms at 164 G §. In the 
sodium experiment, the scattering length a was observed 
to vary dispersively as a function of the magnetic field B, 
in agreement with the theoretical prediction 0]: 
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where ao is the off-resonant scattering length, and A char- 
acterizes the width of the resonance. 

In this Letter we report on the observation of a broad 
Feshbach resonance in sodium in the F = l,mj? = —1 
state at 1195 G, and we investigate the strong inelas- 
tic processes accompanying all three sodium resonances, 
which result in a rapid loss of atoms while approaching or 
crossing the resonances with the external magnetic field. 
The losses show an unpredicted dependence on the ex- 
ternal magnetic field and impose strong constraints on 
future experiments exploiting the tunability of the scat- 



tering length. 

The experimental set-up is very similar to that de- 
scribed in H]. Magnetically trapped condensates in the 
F = l,mj? = —1 state were transferred into an opti- 
cal trap consisting of a focused far off-resonant infrared 
laser beam Q. The atoms could be transferred to the 
mp = +1 state by an rf-pulse in a 1 G bias field. For 
the studies of the Feshbach resonances, bias fields of up 
to 1500 G were applied with ramp speeds of up to 1000 
G/ms. Inhomogenities in the bias field exerted a force 
mainly in the axial direction. To prevent atoms from es- 
caping, green light from an argon-ion laser was focused 
into two light sheets at the ends of the cigar shaped con- 
densates, forming repulsive light barriers. The resulting 
trapping frequencies were around 1500 Hz radially and 
150 Hz axially. The trap losses were studied by ramping 
the bias field with different ramp speeds to various field 
values near the resonances. The atoms were probed in 
ballistic expansion after suddenly switching off the opti- 
cal trap. The magnetic field was switched off 2 ms later 
to ensure that the high field value of the scattering length 
was responsible for the acceleration of the atoms. After 7 
- 25 ms free expansion the atoms were optically pumped 
into the F = 2 state and observed in absorption using the 
cycling transition in a small bias field of 0.5 G. 

The number of atoms N, the scattering length a and 
the mean density (n) could be obtained from the absorp- 
tion images. N is calculated from the integrated optical 
density. The mean field energy 2nh 2 a (n)/m is converted 
into kinetic energy mv^ ms /2 after switching off the trap. 
The root-mean-square velocity v^ ms can be extracted 
from the size of the cloud after the time of flight. As 
discussed in Ref. Q , the mean density (n) is proportional 
to N(Na)~ 3 / 5 for a three-dimensional harmonic oscillator 
potential, and thus the scattering length a scales as 
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Normalized to unity far away from the resonance, v^ ms /N 
is equal to a/ao. The proportionality factor in eqn.(|^) 
involves the mean trapping frequency which was not ac- 
curately measured. However, we can obtain absolute val- 
ues by multiplying the normalized scattering length with 
the theoretical value of a^. The off-resonance scattering 
length ao = 2.75 nm of sodium at zero field increases to 
the triplet scattering length ao — fly = 3.3 nm || at the 
high fields of the resonances. Using this value, the mean 
density is obtained from the root-mean-square velocity, 
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The peak density no is given by no 
potential. 
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Figure 1: Number of atoms N, normalized scattering 
length a/ao & n d mean density (n) versus the magnetic 
field near the 907 G and 1195 G Feshbach resonances. 
The different symbols for the data near the 907 G reso- 
nance correspond to different ramp speeds of the magnetic 
field. All data were extracted from time-of-flight images. 
The errors due to background noise of the images, ther- 
mal atoms, and loading fluctuations of the optical trap 
are indicated by single error bars in each curve. 

The experimental results for the 907 G and 1195 G res- 
onances are shown in Fig. 1. The number of atoms, the 
normalized scattering length and the density are plotted 
versus the magnetic field. The resonances can be identi- 
fied by enhanced losses. The 907 G resonance could be 
approached from higher field values by crossing it first 
with very high ramp speed. This was not possible for the 



1195 G resonance due to strong losses, and also not for 
the 853 G resonance due to its proximity to the 907 G 
resonance and the technical difficulty of suddenly revers- 
ing the magnetic field ramp. For the 907 G resonance the 
dispersive change in a can clearly be identified. The solid 
lines correspond to the predicted shapes with width pa- 
rameters A = 1 G in Fig. (b), and A = — 4 G in Fig. (e). 
The negative width parameter for the 1195 G resonance 
reflects the decreasing scattering length when the reso- 
nance is approached from the low field side. The uncer- 
tainties of the positions arc mainly due to uncertainties 
of the magnetic field calibration. 

The time dependent loss of atoms from the condensate 
can be parameterized as 
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where K{ denotes an z-body loss coefficient, and (n) the 
spatially averaged density. In general, the density n de- 
pends on the number of atoms N in the trap, and the loss 
curve is non-exponential. One-body losses, e.g. due to 
background gas collisions or spontaneous light scattering, 
are negligible under our experimental conditions. An in- 
crease of the dipolar relaxation rate (two-body collisions) 
near Feshbach resonances has been predicted How- 
ever, for sodium in the lowest energy hyperfine state F = 
l,mp = +1 binary inelastic collisions are not possible. 
Collisions involving more than three atoms are not ex- 
pected to contribute. Thus the experimental study of the 
loss processes focuses on the three-body losses around the 
F = 1, mp = +1 resonances, while both two- and three- 
body losses must be considered for the F = l,mp = — 1 
resonance. 

Figs. 1 (c) and (f ) show a decreasing or nearly constant 
density when the resonances were approached. Thus the 
enhanced trap losses can only be explained with increas- 
ing coefficients for the inelastic processes. The quantity 
N/N(n 2 ) is plotted versus the magnetic field in Fig. 2 for 
both the 907 G and the 1195 G resonances. Assuming that 
mainly three-body collisions cause the trap losses, these 
plots correspond to the coefficient K 3 . For the 907 G 
(1195 G) resonance, the off-resonant value of N/N (n 2 ) is 
about 20 (60) times larger than the value for K% measured 
at low fields fij. Close to the resonances, the loss coeffi- 
cient strongly increases both when tuning the scattering 
length larger or smaller than the off-resonant value. Since 
the density is nearly constant near the 1195 G resonance, 
the data can also be interpreted by a two-body coeffi- 
cient K2 = N/N{n) with an off-resonant value of about 
30 x 10~ 15 cm 3 /s, increasing by a factor of more than 50 
near the resonance. The contribution of the two processes 
cannot be distinguished by our data. However, dipo- 
lar losses are much better understood than three-body 
losses. The off-resonant value of K2 is expected to be 
about 10 -15 cm 3 /s in the magnetic field range around 
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the Feshbach resonances, suggesting that the three-body 
collisions are the dominant loss mechanism also for the 
1195 G resonance. 
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a dimer by a third atom as the inverse process of recombi- 
nation, Moerdijk et al. |l2|] suggested that the recombina- 
tion rate is proportional to a 2 . Although the assumptions 
of those theories might not be fulfilled near Feshbach res- 
onances, they raised the hope that loss rates should not 
increase in the region of the Feshbach resonance where the 
scattering length is small. However, for scattering lengths 
smaller than an a substantial increase of the coefficient 
K3 was observed. Thus, these measurements show the 
need for a more accurate theoretical treatment of ultra- 
cold three-body collisions. 

Another way to characterize the losses is a rapid sweep 
across the resonances. We determined the fraction of 
atoms which were lost in sweeping through the 853 G 
resonance and the 907 G resonance at different ramping 
speeds [ pL~3| . For this, the optical trap and the magnetic 
field were suddenly switched off, either before or after 
crossing the resonance, and the number of atoms was de- 
termined from absorption images as before. 



Figure 2: Rate coefficients for inelastic collisions near Fes- 
hbach resonances. N/N{n 2 ) is plotted versus the mag- 
netic field. The time derivative N was calculated from 
neighboring points without smoothing the data. 



The data analysis following eqns.(Q) and (Q) is based 
on two major assumptions: the Thomas-Fermi approx- 
imation in which the kinetic energy is small compared 
to the mean-field energy, and the assumption that the 
atoms maintain the equilibrium density distribution dur- 
ing the change of the scattering length. The Thomas- 
Fermi approximation is well justified for all data points 
with a/an > 1. For a/an < 1 close to the resonances, the 
mean field energy of the condensates is only a few times 
larger than the level spacing in the trap. This leads to 
an overestimate of the density and thus to an underesti- 
mate of K3 for the data points with smallest scattering 
length. The ramp speeds of the magnetic field were cho- 
sen low enough to ensure adiabaticity, characterized by 
the condition a/a <Wj where u>i — 2"KVi are the trap- 
ping frequencies. Only for the data point closest to the 
1195 G resonance (in parentheses) this condition is not 
fulfilled. Indeed, the scattering length a/ao (Fig. lb) and 
the quantity N/N(n 2 ) (Fig. 2a) are independent of the 
ramp speeds, supporting the assumption of adiabaticity. 

The observed increase and the magnetic field depen- 
dence of the three-body collision rate is not accounted for 
by any theory. Two theoretical treatments suggested that 
the rate coefficient K3 should vary monotonically with the 
scattering length following some power law. Fedichev et 
al. jlO| derived the universal relation K3 = 3.9ha 4 /2m. 
This prediction was in fairly good agreement with mea- 
surements in rubidium 1 1 1 and sodium in low magnetic 
fields, even though it is not clear that the assumptions are 
valid in these experiments. By considering the breakup of 
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Figure 3: Fraction of lost atoms after crossing the Fesh- 
bach resonances at 853 G and 907 G for different inverse 
ramp speeds of the magnetic field. 

Fig. 3 shows the fraction of lost atoms as a function 
of the inverse ramp speed. Using the width of 1 G for 
the 907 G resonance, this implies that 70 % of all atoms 
are lost in one microsecond. Across the 853 G resonance 
70 % of all atoms are lost at a ramp speed of 1 G/400^s. 
Assuming a universal behavior near Feshbach resonances, 
this implies that the strength of the 853 G resonance is 
four hundred times weaker than of the 907 G resonance. 
Since the width A of the Feshbach resonance is propor- 
tional to the coupling strength between the two involved 
molecular states §, a width of 0.0025 G for the 853 G 
resonance is estimated. A theoretical prediction for the 
widths of the 907 G and the 853 G resonances are 1 G 
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and 0.01 G, respectively @. 

The loss of 70 % of the atoms in only one microsec- 
ond cannot be explained by the usual picture of inelastic 
collisions. Due to the very small kinetic energy of Bose- 
Einstein-condensed atoms, which is of order ^— (-sr) 2 / 3 
JL5| , the travel distance of an atom in 1 /is is less than 
1 nm under our conditions, much smaller than the mean 
distance between the atoms, d w n -1 / 3 w 100 nm. Loss 
by two- or three-body collisions should then be limited 
to the small fraction of atoms which happens to be very 
close to other atoms. Possible explanations for the ob- 
served large losses are the divergence of the scattering 
length which allows for extremely long range interactions 
between atoms, the formation of a molecular condensate, 
as recently proposed jlq| , or the impulsive excitation of 
solitons |l7j when the region of negative scattering length 
is crossed. All suggested explanations require new con- 
cepts in many-body theory. 

In conclusion, we have reported on the observation of 
a new, F = l,mj? = —1 Feshbach resonance at 1195 G 
in sodium, and we have investigated the enhanced trap 
losses accompanying all three resonances observed so far. 
The scattering length could be altered in the range of 
0.2 < a/a,Q < 5. The three-body loss coefficient K$ in- 
creased for both larger and smaller scattering length by 
up to more than three orders of magnitude. Those losses 
show that the density of the sodium samples must be re- 
duced well below 10 14 cm -3 for further studies exploiting 
the tunability of the scattering length near Feshbach res- 
onances. 

Our sweep experiments revealed that approaching a 
Feshbach resonance from the high magnetic field side is 
strongly affected by trap loss while crossing the resonance. 
The new rriF = —l resonance at 1195 G is well suited for 
studies of Bose-Einstein condensates with zero or negative 
scattering length, since this region is on the low field side 
of the resonance and can be directly approached without 
crossing any resonance. 

So far, the physics of gaseous Bose-Einstein conden- 
sates has been very well described including only binary 
interactions between the atoms |l5j . Feshbach resonances 
might open the possibility to study physics beyond this 
approximation which breaks down when the scattering 
length diverges, and also for very small scattering lengths 
when one has to consider higher order terms in the atomic 
interactions. Feshbach resonances may also lead to new 
inelastic processes. The observed losses of atoms near 
Feshbach resonances indicate molecular and many-body 
physics which is not yet accounted for by any theory. 
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